Glass-Like Heat Conduction in High-Mobility Crystalline Semiconductors 
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The thermal conductivity of poly crystalline semiconductors with type-I clathrate hydrate crystal 
structure is reported. Ge clathrates (doped with Sr and/or Eu) exhibit lattice thermal conductiv- 
ities typical of amorphous materials. Remarkably, this behavior occurs in spite of the well-defined 
crystalline structure and relatively high electron mobility (~ 100 cm 2 /Vs). The dynamics of dopant 
ions and their interaction with the polyhedral cages of the structure are a likely source of the strong 
phonon scattering. 



PACS numbers: 66.70.+f, 63.20.-e, 63.20.Pw, 63.50,+x 



The low-temperature thermal conductivity (k) of 
glasses has universal features that have been the subject 
of experimental and theoretical investigation for many 
years Q; at the lowest temperatures (T < 1 K), k oc T 2 , 
and for 4 K< T <20 K, k is nearly temperature inde- 
pendent (the "plateau"). The T 2 behavior is attributed 
to phonon scattering by a broad distribution of two-level 
systems, whereas the plateau arises from a strong de- 
crease in the phonon mean free path with increasing fre- 
quency, e.g. due to Rayleigh scattering from variations 
in density or bond length. Though typical of amorphous 
materials, similar thermal conductivity behavior, with k 
approaching theoretical minimum values [0 , has been ob- 
served in some disordered crystalline dielectric materials 
HQ ■ A common characteristic of these latter systems is 
the presence of atoms or molecular groups having two or 
more semi-stable positions, and the absence of long-range 
correlation between their positions or orientations. Res- 
onant scattering of phonons by these "rattling" entities 
can also produce a plateau in k || . 

The search for high-mobility semiconductors with 
glass- like thermal properties -"phonon glass, electron sin- 
gle crystals" (PGEC's) - was initiated by Slack M, and 
motivated by thermoelectric applications @ which could 
be substantially improved by the higher figure of merit 
such hypothetical materials would possess. It was pos- 
tulated that the group IV elements forming clathrate- 
structured compounds were potential PGEC's M. In 
this Letter we report the realization of PGEC behav- 
ior in crystalline Ge clathrates. Trends in the behav- 
ior of the low-temperature lattice thermal conductivity 
for Si, Ge and Sn clathrates indicate that the glass- 
like heat conductivity arises from resonant scattering of 
phonons by vibrations of the dopant ions in their poly- 
hedral cages. The relatively high charge-carrier mobility 
of these clathrates indicates that the dopant-cage inter- 
actions do not substantially degrade the electronic prop- 
erties. 



The sample preparation has been described previ- 
ously 0. The specimens are polycrystalline Si and 
Ge clathrates having the type-I hydrate crystal struc- 
ture [fHH , doped with divalent barium (Ba 2+ ), strontium 
(Sr 2+ ), and europium (Eu 2+ ) ions, and Sn clathrates of 
the same structure, doped with monovalent Cs + . The 
cubic crystal lattice is a tetrahedral network of Si, Ge or 
Sn which forms periodic voids or " cages" of 20- and 24- 
coordinated polyhedra in a 1:3 ratio, respectively. The 
dopant ions enter the lattice interstitially, residing inside 
these large cages. X-ray diffraction (XRD) spectra from 
powdered clathrate samples revealed sharp lines with 
no phases other than that of type-I clathrate-hydrate 
(Fm3n). Trace amounts of Si were detected in the Si- 
clathrate sample. Electron beam microprobe analysis of 
a polished surface of each sample confirmed the x-ray 
results and yielded stoichiometric percentages for each 
sample. In the case of the Ge-clathrates, scanning tunnel- 
ing microscopy showed a defect-free structure within the 
grains. Table I lists some of the physical and electronic 
properties of these phase-pure, dense polycrystalline sam- 
ples. Gallium, which substitutes for Si or Ge, was added 
to compensate for the charge of the divalent dopants. 
Charge neutrality in the CssSn44 compound implies two 
Sn vacancies per unit cell. A preliminary report Q on 
specimens with nominal composition SrgGaigGe3o indi- 
cates that the carrier concentration and mobility vary 
with composition. Their relatively high mobilities (up to 
100 cm 2 /Vs) in comparison to other thermoelectrics JlC[ j 
suggest the possibility that the Ga may be ordered in 
these materials, an issue currently under investigation. 

The samples were cut with a wire saw into paral- 
lelepipeds with approximate dimensions 5x2x2 mm 3 . 
Steady-state thermal conductivity measurements at T > 
4 K were conducted in a radiation-shielded cryostat 
using a single heater and 25/im-diameter differential 
chromel-constantan thermocouple. The electrical resis- 
tivity and thermoelectric power were measured simul- 
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taneously; room-temperature values are listed in Table 
I. Heat losses due to radiation and conduction through 
the leads were measured in separate experiments and the 
data corrected. These losses were 10-15% near room tem- 
perature and < 5% for T <120 K. The absolute accuracy 
in k is ±10%, limited by uncertainty in the thermocouple 
junction separation. Measurements at T < 4 K were con- 
ducted in a dilution refrigerator using a similar technique 
with two thermometers. 

Figure 1 shows lattice thermal conductivities, 
K e , in the range 5 K< T <300 K for the five clathrates 
listed in Table I. The electronic thermal conductivity, 
K e , was estimated using the electrical resistivities and 
Wiedemann-Franz law with ideal Lorenz number, Lq = 
2.44 x 10~ 8 WS1/K 2 . Also shown is k for amorphous ger- 
manium (a-Ge) and a-Si0 2 @. The X=295 K values of 
Kg for all of the specimens are small, lying between that 
of a-SiC>2 and a-Ge, and approaching K m i n for Ge, the 
minimum theoretical value calculated using the model of 
Slack Q| , and taking the minimum mean free path (mfp) 
of the acoustic phonons as one-half their wavelength jj) . 
The Sn clathrate exhibits a temperature dependence typ- 
ical of crystalline insulators, with K g increasing with de- 
creasing T approximately as 1/T, the signature of prop- 
agating phonons scattered by anharmonic interactions. 
For the Si-clathrate, n g also increases with decreasing T, 
but not as strongly, indicating additional scattering at 
low temperatures. The Ge clathrates have n g suppressed 
below that of the Sn clathrate by more than an order of 
magnitude at low T. Most striking are the plateaus or 
minima in n g (T) for the Ge clathrates, similar to that 
of a-Si02 and other amorphous materials in the range 
4 K< T <30 K. This similarity extends to T < 1 K as 
shown in Fig. 2 for the SrgGaigGe3o specimen; n g oc T 2 
is observed down to T—60 mK. Our previous estimates of 
grain-boundary scattering for the Ge clathrates M and 
preliminary k measurements on Ge clathrates with sim- 
ilar composition and mm-sized grains fio)] , indicate that 
grain boundaries are not the source of the glass-like k in 
these materials. 

To interpret the n g data in the absence of experimen- 
tal studies of the vibrational spectra for these semicon- 
ducting clathrates, it is instructive to consider prior work 
on the thermal conductivity Jl2[ and vibrational prop- 
erties |l3) of type I clathrate hydrates (ice clathrates). 
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FIG. 1. Lattice thermal conductivity vs. temperature for 
the five clathrate specimens listed in Table I. The dashed and 
dotted curves are for a-Sf0 2 and a-Ge (Ref. @), 

and the 

solid curve is the calculated minimum thermal conductivity 
of Ge (see text). 

comprised of water molecules, which encapsulate gas 
atoms or molecules (guests). Their thermal conductiv- 
ities are anomalously low and glass like. The qualitative 
features of the guest-host interactions for these materi- 
als are expected to be similar to those of semiconduct- 
ing clathrates. Neutron scattering studies and molecu- 
lar/lattice dynamical calculations |Q have established 
that localized guest vibrations interact strongly with the 
host acoustic phonon branches in the clathrate hydrates. 
This interaction arises as an avoided crossing of dispers- 
ing phonon branches and non-dispersing guest vibrations 
of the same symmetry, and yields a resonant damping of 
heat-carrying phonons at certain points in the Brillouin 
zone. The guest translational vibration frequencies in- 
crease as the size difference between guest and host cage 
decreases, a natural consequence of stronger restoring 
forces for the guests. Generally, one expects two promi- 
nant resonances, with the lower (higher) frequency corre- 
sponding to guest vibrations in the larger (smaller) host 
cages. 



TABLE I. The five compounds of this study with atomic percentages from electron-beam microprobe analysis, lattice 
parameters, ao, average grain sizes, room-temperature electrical resistivities, p, and absolute Seebeck coefficients, S. 



compound 


Elemental atomic % 


o (A) 


grain size (/im) 


p(mSlcm) 


S(/xV/K) 


CS8Sn44 


15.6Cs/84.4Sn 


12.113 


13.0 


59.4 


-304 


BasGaigSiso 


15.5Ba/27.5Ga/57.0Si 


10.554 


9.6 


2.04 


-66 


Sr g Gai6Ge3o 


14.8Sr/30.2Ga/55.3Ge 


10.731 


9.6 


12.8 


-313 


EusGai6Ge3o 


14.2Eu/28.3Ga/57.6Ge 


10.711 


5.4 


2.52 


-152 


Sr 4 Eu4Gai6Ge3o 


8.2Sr/6.9Eu/28.6Ga/56.4Ge 


10.726 


12.2 


1.00 


-88 
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FIG. 2. Lattice thermal conductiviy measured from 60 mK 
to 100 K for SrsGai6Ge3o with a fit (solid curve) to the model 
discussed in the text (see also Fig. 3). The dashed line is for 
a-Si02, and the solid curve labelled Kmin is the calculated 
minimum thermal conductivity of Ge. 

Consider the relevant guest and host dimensions for 
our specimens. The radii of the Ba 2+ , Eu 2+ and Sr 2+ 
ions are 2.07, 1.91 and 1.90 A, respectively, from x- 
ray crystallographic data of the respective tclluridc com- 
pounds with the NaCl-type structure jl5| . These dimen- 
sions may not be accurate for the Si or Ge-clathrates 
where the ions reside in large polyhedral "cages" , but 
the trend in radii should be similar. The S120 and Si24 
cage radii are estimated to be 2.05 and 2.23 A, respec- 
tively, from XRD data on Si-clathrates doped with alkali 
metals jl6|. XRD also indicates no significant influence 
of Ga substitution on the average size of the cages. The 
Ge2o and Ge24 cages have radii of 2.15 and 2.32 A, respec- 
tively Jig] . Thus we expect that Eu 2+ and Sr 2+ in the 
Ge compounds are "looser" in their cages than are Ba 2+ 
ions in the Si compound. We infer that the ionic radius 
for Cs + in the Sn clathrate is approximately the same 
as the cage size from the following analysis. Using data 
for isostructural K- and Rb-doped Sn clathrates from the 
literature p"^ ] we plot lattice constant vs alkali-metal ion 
radius [Q, and observe that the measured lattice con- 
stant for CsgSn44 is larger than that implied by linear 
extrapolation on this plot. This suggests that Cs + in the 
cages of CssSn44 causes a lattice expansion. 

From these observations it is apparent that the mag- 
nitude of the T = 10 K lattice thermal resisitivity, 
1/k 9 , correlates with the relative guest/host-cage size 
mismatch. This motivates our proposal that the plateaus 
or dips in K g for the Ge clathrates are associated with 
resonant scattering via dopant vibrations. The ab- 
sence of this behavior in the Sn clathrate and weaker 
effect in the Si clathrate suggest that the vibrational 
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FIG. 3. Difference in the lattice thermal resistivity of each 
Ge clathrate from that of the Sn clathrate. The frequency 
scale (upper ordinate) is for dominant thermal phonons 
(hu — 3.8fcflT). The inset shows the inverse scattering length 
at T = 10K for the resonant terms (see text) used to fit the 
data for SrsGai6Ge3o (solid curve in Fig. 2). The parameters 
were: G x = 3.8 x lOV^s^K" 2 , C 2 = 1.7 x lOV^s^K- 2 , 
o>i = 1.7THz, u 2 = 6.3THz, and 71 = 72 = 0.8. 



frequencies of Cs + and Ba 2+ , respectively, are sufficiently 
high in energy that a strong interaction with the acous- 
tic phonon branches of the host does not occur for these 
compounds. To further explore the resonance scatter- 
ing within this scenario we use the Sn-clathrate data as 
a reference, computing the difference in thermal resis- 
tivity for each of the other clathrates from that of the 
Sn clathrate, AW — l/n g — l//%sn- These results are 
plotted in Fig. || and, if this interpretation is correct, rep- 
resent a measure of the thermal resistivity due to reso- 
nance scattering in the Si and Ge clathrates. We see that 
AW for the Ge clathrates is peaked in the temperature 
range 10-20 K. Within the dominant phonon approxima- 
tion {Tiujdom — 3.8ksT for constant mfp) this corresponds 
to phonon frequencies oJdom ~ 5 — lOTHz (upper ordinate 
in Fig. 3). The asymmetry of the maxima for the Sr- and 
Eu-doped specimens suggests the presence of a second, 
smaller resonance peak at lower frequency. 

The K g oc T 2 behavior at low temperatures implies a 
phonon mfp inversely proportional to frequency, l(u>) oc 
lo^ 1 . This follows from the kinetic theory expression, 
K g = Cvl/3 (C and v are the lattice specific heat 
and sound velocity, respectively), a Debye specific heat, 
C oc T 3 , and to oc T for thermal phonons. The localized 
tunnel systems (TS) that give rise to this phonon scat- 
tering in amorphous materials Q also yield anomalous 
features in the specific heat and elastic properties in the 
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same temperature range. It is plausible that TS's are 
the source of the scattering in the clathrates. For exam- 
ple, they could correspond to different positions of the 
dopant ions in their cages; a broadened density of these 
states might be induced by the random distribution of 
the dopants. Measurements of specific heat and elastic 
properties are needed to address this issue. 

At T < 40 K, where n g 3> n m i n and phonon trans- 
port is an appropriate description, the behavior of K g for 
the Ge specimens can be simulated by integrating the 
kinetic theory expression, n g = (u/3) J C{u>)l(ui), with 
a Debye specific heat and a mfp that is a sum of terms 
representing TS, resonant, and Rayleigh (R) s cattering, 

l(u) = (Irl + l^es + IrY 1 + Imin, with |§§(||l), 

i-l = A(haj/k B )t&rih{huj/2k B T) 

+ (A/2)(k B /hu; + B- 1 T- 3 y 1 , 

l~ e \ = £ C^l^l c 2 ) 2 + >y iU * U % 
i 

I], 1 = Dh 4 uj 4 /k B . 

The constants A and B are related to microscopic vari- 
ables describing the TS model |pl| . The phenomeno- 
logical resonance terms are of the form employed pre- 
viously to describe phonon scattering in ionic crystals 
p2fl . The lower limit on I is assumed to be a constant, 
Imin- Using v = 3370 m/s and B — 360 K appropriate 
for diamond-structured Ge ^3|, the solid curve in Fig. 2 
was generated for Sr 8 Gai 6 Ge3o using A = 1.38 x 10 4 
m -1 K _1 , B = 1.5 X 10~ 3 K" 2 , l min = 3 A, D = 1 
m _1 K -4 and a sum of two resonance terms shown in the 
inset of Fig. 3. The TS parameters are comparable to 
those found for many amorphous solids | pl[ . The l~ e \ 
employed is not unique and the values of to\ and 102 are 
sensitive to the value of D in the range characteristic 
of crystalline materials, < D < 5 m _1 K -4 (two or- 
ders of magnitude smaller than in glasses |fl]]). Within 
this range of D, we achieved good fits for resonance fre- 
quencies, 0.6 THz< lux < 2 THz and 3 THz< u 2 < 20 
THz. The data for the other Ge specimens could be 
fitted using the same TS parameters by simply adjust- 
ing the magnitudes of the two resonance terms without 
substantially altering the resonance frequencies. Thus a 
sum of two resonances has the minimum complexity re- 
quired to reproduce the range of n g behavior we have 
observed. We expect a "freeze-out" of the guest-ion rat- 
tling at low temperatures; this should give rise to static 
disorder in bond lengths associated with a distribution of 
the guest-ion positions within their cages. Small energies 
(< 1 K) presumably characterize the separation between 
these states, and are the origin of the TS. 

These unique materials possess interesting and poten- 
tially useful properties that call for further investigation. 
As one of only a few crystalline materials having glass- like 
thermal conduction, other thermal and elastic properties 



at T < 1 K are of particular interest. 
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